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The electromagnetic form factors of the three-nucleon bound states were calculated in Complete 
Impulse Approximation in the framework of the Covariant Spectator Theory for the new high- 
precision two-nucleon interaction models WJC-1 and WJC-2. The calculations use an approximation 
for the three-nucleon vertex functions with two nucleons off mass shell. The form factors with WJC-2 
are close to the ones obtained with the older model W16 and to nonrelativistic potential calculations 
with lowest-order relativistic corrections, while the form factors with the most precise two-nucleon 
model WJC-1 exhibit larger differences. These results can be understood when the effect of the 
different types of pion-nucleon coupling used in the various models is examined. 
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I. INTRODUCTION 

The electromagnetic form factors of nuclei provide im- 
portant information about their internal structure. They 
have been used extensively in order to test models of the 
nuclear dynamics and of the associated electromagnetic 
currents. As electron scattering experiments, such as the 
ones performed at Jefferson Lab, reach larger and larger 
values of the momentum transferred by a virtual photon 
to the struck nucleus, it becomes increasingly important 
to incorporate the requirements of special relativity in a 
reliable way into the theoretical description of the pro- 
cess. 

The Covariant Spectator Theory (CST) Q was de- 
signed as a manifestly covariant theory, especially suited 
for the description of few-nucleon problems. In a re- 
cent paper Q, we presented the first CST calculations 
of the electromagnetic three-nucleon (3N) form factors 
in Complete Impulse Approximation (CIA), which is de- 
fined as the complete CST 3N current [|[ except for inter- 
action currents, i.e., diagrams where the photon couples 
to an intermediate interacting two-nucleon (NN) system. 
However, the term "impulse approximation" can be mis- 
leading because it depends on the framework used. 

For instance, in a very successful approach used by the 
Pisa- Jlab collaboration 0, HI , the dynamics is based on 
the nonrelativistic Schrodinger equation, and relativistic 
corrections are added perturbatively. We call the corre- 
sponding impulse approximation with relativistic correc- 
tions "IARC" . In this framework, two- and three-body 
interaction currents are later added to the IARC re- 
sults. These include first-order 'jttNN contact interac- 
tions which are equivalent to the already at the CIA level 
automatically — and to all orders — included "Z-graphs" 
in the CST. This is an example of a more general observa- 
tion: what counts as interaction current in one approach 



may be part of the impulse approximation in another. 

In Ref. @, our focus was to study the model depen- 
dence of the electromagnetic form factors in CST. We 
performed calculations for a family of closely related rel- 
ativistic two-nucleon interaction models and found that 
the CST results behave very reasonably. In most cases, 
a direct comparison of our CIA results with experimen- 
tal data is not useful because we expect interaction cur- 
rents to be significant. However, the comparison with 
IARC results is instructive, and it appears that the sur- 
prisingly close agreement between the two approaches 
(at least for Q < 4 fm _1 ), when models are compared 
that yield the same 3N binding energy, is no coinci- 
dence. The main reason seems to be that all CST mod- 
els used in the comparison employ pseudovector cou- 
pling for the pion-nucleon vertices. This kind of coupling 
suppresses negative-energy states (corresponding to Z- 
graphs), which are included in CIA but not in IARC. 
It is therefore understandable that no large differences 
between the two calculations emerge, as long as other 
aspects of the dynamics in the two approaches are com- 
parable. 

It would be interesting to submit this interpretation to 
a test. One only needs to perform two calculations in CIA 
with two NN models that are as similar as possible in 
their ability to describe the NN data and the 3N binding 
energy. One of them should be based on pure pseudovec- 
tor pion-nucleon coupling, while the other should include 
an admixture of pseudoscalar pion-nucleon coupling and 
thus increase the weight of Z-graphs. If the above inter- 
pretation is correct, the model with pure pseudovector 
coupling will be close to the IARC result, while there 
should be larger deviations in the case of the model with 
some pseudoscalar coupling. 

We are indeed in a position to perform this test. In 
a recent paper Q, we published two realistic CST mod- 
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els for the neutron-proton interaction, both of which de- 
scribe the np scattering observables with xV-^data ~ 1 
for the most recent 2007 data base. The first model, 
WJC-1, based on the exchange of 8 bosons and fitted with 
27 adjustable parameters, features a mixture of pseu- 
dovector and pseudoscalar pion-nucleon coupling. The 
second model, WJC-2, based on the exchange of 6 mesons 
and with only 15 adjustable parameters, uses pure pseu- 
dovector pion-nucleon coupling. The two models can be 
considered to be essentially on-shell equivalent, and both 
reproduce also the experimental value of the triton bind- 
ing energy of 8.48 MeV. 

There is, however, one obstacle to performing the CIA 
calculations with models WJC-1 and WJC-2: Some of 
the diagrams that comprise the CIA 37V current depend 
on the 3iV vertex function with two nucleons off mass 
shell. A computer code for the calculation of these vertex 
functions for the new models WJC-1 and WJC-2 is at 
present in development, but not yet ready to be used in 
the calculation of the 3N form factors. 

This obstacle can be overcome if we apply an approx- 
imation in which 3N vertex functions with two nucleons 
off mass shell are appropriately replaced by vertex func- 
tions with only one nucleon off mass shell. We can test 
the quality of this approximation, which we call "CIA- 
0" , by applying it to one of the models previously used 
in Q and comparing the approximate form factors to the 
respective full CIA result. 

Of course, if it turns out that CIA-0 is a reliable ap- 
proximation, the 37V form factors obtained from the re- 



alistic models WJC-1 and WJC-2 will be of high interest 
by themselves, not just as a means to evaluate the effect 
of Z-graphs. While the family of models used in 0] gives 
a good description of the NN data, they cannot compete 
in precision with our new models. 

This paper is divided into four sections. After the in- 
troduction, Section [II] describes the 37V current and de- 
fines the CIA-0 approximation. Section HTT1 presents and 
discusses the numerical results obtained, and in Section 
IIVI we draw our conclusions. 



II. THE 3N CURRENT IN CIA AND CIA-0 

The complete form of the electromagnetic 3N current 
in CST was derived in Ref. Q, and used in Q for the 
first time to calculate the 3iV form factors in CIA. Fig- 
ure [T] displays the complete current, and CIA is defined 
through diagrams A to F. 

We denote the photon four-momentum by q, and we 
label the nucleon four-momenta ki such that always 
k\ = k 2 = to 2 , where to is the nucleon mass. For the 
cases where a nucleon absorbs a photon, we introduce 
the notation kf~ = fcj ± q. The momentum k% is not an 
independent variable; in the CST, the energy-momentum 
four-vector is conserved and it is determined through the 
momenta of nucleons 1 and 2 and the total 3N momenta 
P t in the initial and P' t = Pt + q in the final state. 

The 3iV current in CIA is given in algebraic form by 
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where E(k) = V to 2 + k 2 , V12 is a permutation operator 
which interchanges particles 1 and 2, £ is a phase with 
( = +1(— 1) for bosons (fermions), Gp'p(k) is the propa- 
gator of an off-shell nucleon with four-momentum k, and 
ja'a(k',k) is the single nucleon current for off-shell nu- 
cleons with incoming (outgoing) four-momentum k (k 1 ). 
Summation over repeated Dirac indices is implied. The 
3N vertex functions T arc solutions of Faddeev-type CST 
integral equations [7J] and were obtained numerically for 
the NN interaction models considered here [6, 8]. The 
"relativistic wave functions" are defined as 

*A 1 A 2 a(fci,fc2;-Pt) = G aa '(k 3 )Tx 1 x 2a '{h,k 2 ;Pt) , (2) 
and we use a shorthand for the contraction of Dirac in- 



dices with nucleon helicity spinors with hclicity A.;, 



^Xi\ 2 a'(ki, k 2 ; P t ) = 

u ai (ki, Ai)w Q2 (fc 2 , A 2 )r (k u k 2 ;P t ). (3) 



In the second and third line of Eq. |T]), corresponding 
to diagrams (B+C+E+F) of Fig. [TJ the vertex function 
appears with two nucleon momenta off mass shell. The 
solutions of the CST equation for the 3N bound state 
have only one nucleon (nucleon 3, by convention) off mass 
shell, but one can obtain vertex functions with two off- 
shell particles through an iteration of the 3iV equation 
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FIG. 1: (Color online) The electromagnetic 3iV current in CST for elastic electron scattering from the 3iV bound state. A cross 
on a nucleon line indicates that the particle is on mass shell. Diagrams (A) to (F) define the complete impulse approximation 
(CIA), in which the photon couples to single nucleons, which can be off-shell (A and D) or on-shell before or after the photon- 
nucleon vertex (B, C, E, and F). The approximation denoted CIA-0 replaces the vertex function with two nucleons off mass 
shell in diagrams (B), (C), (E), and (F) by a vertex function with only one nucleon off mass shell. The interaction diagrams 
(G) to (J) describe processes in which the photon couples to two-body currents associated with the two-nucleon kernel. 



with an off-shell two-nucleon scattering amplitude, 

m d 3 ki, 
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Here, Mp a y^ a i (k 2 , k' 2 ; P23) is the scattering amplitude 
of nucleons 2 and 3 with total pair momentum P23 > and 
the final state momentum of nucleon 2, k 2 , is off mass 
shell, while its initial-state momentum, k 2 , is on mass 
shell (nucleon 3 is off shell in either state). 

As pointed out in the Introduction, the 37V vertex func- 
tions with both nucleons off-shell in the final state are not 
available at this time for the new NN interactions WJC- 
I and WJC-2. Moreover, it is rather awkward to cal- 
culate and manipulate these double-off-shell vertex func- 
tions numerically, because with one additional continuous 
variable (the off-shell energy of nucleon 2) they occupy 
much more computer storage space and slow down the 
calculations. 

For these practical reasons, we introduce here a simple 
approximation which replaces the vertex functions with 
two nucleons off-mass-shell by others with only one nu- 
cleon off mass shell. 

In order to motivate this approximation, consider for 
instance diagrams (C) and (F) of Fig. [I] The vertex func- 
tion in the initial state, T\ 10a (ki, k 2 ; Pt), depends on the 
off-shell momentum k 2 = k 2 — q. Since k 2 is on mass 
shell, for small photon momenta q the momentum k 2 is 
also almost on mass shell. We may therefore expand the 
vertex function in the off-shell energy of nucleon 2 around 
its on-shell value. If we keep only the zeroth-order term of 



the expansion and eliminate the corresponding negative- 
energy channel of nucleon 2 (with negative p-spin), we 
obtain a known vertex function with two nucleons on 
mass shell. We call this approximation "CIA-0," refer- 
ring to the zeroth-order expansion involved. 

While this approximation is easy to apply, its for- 
mulation is somewhat awkward because of its frame- 
dependence. In our numerical calculations, the 3iV vertex 
function is expressed in terms of variables for nucleons 2 
and 3 which are defined in the rest frame of the (23) pair 
where the CST equation for the two-nucleon scattering 
amplitudes is solved numerically. We can write 

r Xl a {ki,k 2 ;P t ) = r Xl i3a(ki, L(k 23 )k 2 ; P t ) , (5) 

where the Lorentz transformation L{k 23 ) takes the sys- 
tem of nucleons 2 and 3 from its rest frame, where their 
momenta are k 2 and k%, to the 3iV rest frame, where 
their momenta are k 2 = L(k 23 )k 2 and ^3 = L(k 23 )k^ 
and where their total two-body momentum is k 23 — 
k 2 + ks = Pt — fci . 

We define now the four-momentum f 2 to have the 
same three-vector part as k 2 but to be on mass shell, 



(E(k 2 



and we replace the momentum 



k 2 by f 2 in the vertex function (|5j). 

In order to eliminate the negative energy states of nu- 
cleon 2, we first write the propagator of nucleon 2 in 
terms of its form in the pair rest frame, 



Gff'pih ) = Sp lf}l {L{k 23 )) 
I m + f 2 



(k 2 y 



S^ (L(k 23 )) , (6) 
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where is the Dirac space representation of the 

Lorentz transformation L(fc^). 

Now we keep only the component with positive p-spin 
in the pair rest frame, 



m + $ 2 



rn 



A+(f 2 ") 



(7) 



with the positive-energy projector 

m + f 2 



A+(r7) 



2m 



(8) 



The approximation CIA-0 can then be defined as the re- 
placement 



Gpip{k 2 )T Xl p a (ki,k 2 ;P t 
%MM fe 23)) 



E(k 2 -)k 2 - -E(k 2 -)-ie 

xr Al ^(A; 1 , J L(fc 23 )f 2 -;P t ) (9) 

in Eq. (TIj), as well as an analogous replacement for 
^\ 1 0'a(ki,k2;P[)Gf3'f}(k2), which occurs in diagrams 
(B) and (E). 

Note that the projector A+ eliminates negative-energy 
states of nucleon 2 in the two-body rest frame, but this 
does not eliminate all Z-graph contributions from the 
calculation. They are still present through the negative- 
energy states of nucleon 3, and they are also re-generated 
to some extent when the state of nucleon 2 is boosted to 
other frames. 

The approximation ([9]) may look complicated, but it 
is actually easy to implement in our numerical calcula- 
tions. For instance, in the case of diagram C it merely 
amounts to replacing in Eq. (B64) of Ref. [1] the off- 
shell energy p° of nucleon 2 by the corresponding on- 
shell value E(p) in the argument of the partial wave ver- 
tex function C(qp°pM jm\i\2\^p2P^TT z ) 1 and restrict- 
ing the summation over the p-spins of nucleon 2 to the 
positive-energy value p 2 = + only. 

The electromagnetic current for an off-shell nucleon 
can be written in the form 



f N {k',k) =h(k'\k 2 ) F 1N {Q Z ) 1 



+fl ) (k u ,k')F 2N (Q') 



2,i<y^q v 
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+g n (k' 2 ,k 2 )F 3N (Q 2 )A_(k')^A^(k) , (10) 

where fo, fo, and go are nucleon off-shell form factors 
associated with the boson-nucleon vertices, and F±n and 
f 2 jy are the usual electromagnetic Dirac and Pauli form 
factors. Since A_ projects onto negative energy states, 
the form factor F^n belongs to a term that contributes 
only if the nucleon is in a negative-energy state before 
and after the photon-nucleon vertex. We adopt the usual 
convention Q 2 = —q 2 . 



TABLE I: Parameters Fo and k (in fm -1 ) of the scaling func- 
tions F S (Q) of Eq. (|12[) by which the electromagnetic 3N form 
factors are divided in the figures with linear scale. 





charge f.f. 


magnetic f.f. 


Form factor 


Fo 


k 


Fo k 


3 H 


1 


0.760488 


1 0.871664 


3 He 


1 


0.799411 


1 0.912562 


Isoscalar 


1.5 


0.778026 


0.423 0.765425 


Isovector 


0.5 


0.842695 


2.13 0.889873 



The isospin dependence of the electromagnetic form 
factors is, for i = {1,2,3} and the nucleon isospin pro- 



jection T 



F iN (Q 2 ) = F ip (Q 2 )- 



F m {Q 2 ) 



1 



(11) 
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In previous calculations @ , we found that the 3N form 
factors are quite insensitive to the inclusion and varia- 
tions of the off-shell nucleon form factors. Therefore we 
employ in the calculations of this work the simpler on- 
shell nucleon current, with /o = /o = 1 and go = 0. 
For the Dirac and Pauli form factors, we chose the pa- 
rameterization of Galster Q, in order to compare with 
IARC results provided to us by Marcucci [10( who used 
the same parameterization. 

With the CIA-0 approximation in place, the electro- 
magnetic 3Af form factors are calculated numerically 
from the 3N vertex functions, which were obtained by 
solving the 3Af CST equation in helicity partial wave 
form. The applied techniques are described in detail in 
Ref. [|. 



III. PRESENTATION AND DISCUSSION OF 
THE RESULTS 

We calculated the electromagnetic 3A^ form factors for 
three NN interaction models, W16, WJC-1 and WJC-2, 
for momentum transfer up to Q = 9 fm . The results 
are displayed in Figs. and [3J 

Since the form factors fall several orders of magnitude, 
and the traditional log-plots tend to obscure differences 
in some places and overemphasize them in others, we 
divide them by simple scaling functions of the form 



F s (Q)=F e 



-Q/k 



(12) 



Table [J shows the parameters of the scaling function for 
each case. We also list the magnetic moments in Table 
Hi] and the charge and magnetic radii in Tabic [TTTl 

First, we start with a comparison of the curves for 
W16 in CIA and in CIA-0, which clearly demonstrates 
the high quality of the approximation. The differences 
between the exact calculation and the approximation are 
hardly noticeable up to values of Q around 7 fm -1 , and 
in general appear to be insignificant. We may therefore 
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4 5 
Q (fm" 1 ) 

FIG. 2: (Color online) Charge form factors of the 3iV bound states, 3 H (first row), 3 He (second row), and the isoscalar (third 
row) and isovector (fourth row) combinations. In each case, the figure on the left shows the form factor in the traditional 
semi-log plot, while the figure on the right shows the same form factor divided by a scaling function of Eq. (|12[l 2] on a linear 
scale. The solid line is the result for NN model W16 in CIA, the dotted line is the approximation CIA-0 for the same model. 
The dashed line is model WJC-1, and the dash-dotted line is model WJC-2, both in CIA-0. For comparison, the solid line with 
theoretical error bars is the result of an IARC calculation by Marcucci [l(| based on the AV18/UIX potential. All calculations 
employ the on-shell single- nucleon current, with the Galster parameterization of the nucleon form factors The full circles 
represent the experimental data [3, El GJ, GJ, EE EE E3, EE O • 
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TABLE II: Magnetic moments in nuclear magnetons (n.m.). 
The experimental values are from Ref. [201 ] . 

Model M( 3 H) M( 3 He) [is W 

W16 (CIA) 2.544 -1.747 0.400 -2.144 

W16 (CIA-0) 2.543 -1.743 0.400 -2.143 

WJC-1 (CIA-0) 2.441 -1.648 0.396 -2.044 

WJC-2 (CIA-0) 2.525 -1.742 0.391 -2.134 

IARC 2.572 -1.763 0.404 -2.168 

Experiment 2.979 -2.128 0.426 -2.553 



assume that the results for WJC-1 and WJC-2 obtained 
here only in CIA-0 should also be very close to the exact 
CIA result. 

Note that there are caveats to this conclusion: the 
quality of CIA-0 compared to CIA was really tested only 
for W16, a model with a very smooth choice for the 
definition of the kernel (and hence the vertex function) 
when both nucleons off-shell. The WJC models have a 
more complex off-shell structure (corresponding to the 
prescription C discussed in Ref. [(|) and their off-shell ex- 
trapolations will not be as smooth. In addition, WJC-1 
has a mixed pseudoscalar-pseudovector pion-nucleon cou- 
pling, and it is conceivable that the pseudoscalar part of 
this coupling might introduce further differences between 
CIA and CIA-0 to which W16 is not sensitive. Our con- 
clusions must therefore be taken with these particular 
grains of salt. In any case, CIA-0 should be a better 
approximation to CIA at smaller Q, simply because the 
nucleon involved is taken less far off mass shell. 

We turn now to a comparison of the form factors for 
different models of the NN interaction. The figures show 
that the WJC-2 form factors stay close to the ones of 
W16, while, in most cases, WJC-1 begins to deviate 
somewhat already at smaller values of Q. WJC-2 and 
W16 are also close to the IARC results, typically up 
to about Q = 6 fin -1 . This supports the conjecture 
made in the Introduction, namely that the suppression 
of Z-graphs through the use of pseudovector pion-nucleon 
coupling is mainly responsible for the close agreement be- 
tween the CST and IARC. 

Apart from the issue of the type of pion-nucleon cou- 
pling, the CST models include other boson-exchanges 
with off-shell coupling. Most notably, those due to scalar 
isoscalar (ao) and isovector (oi) exchanges have been 
found to have a very strong influence on the quality of 
the NN fits and on the triton binding One might 
expect them to have a strong influence on the 3iV form 
factors as well. 

The results indicate that this is only indirectly the case, 
namely through their effect on the binding energy. When 
the scalar off-shell coupling strength is varied without 
constraining the triton binding energy, the 3N form fac- 
tors show substantial variations (2j. On the other hand, 
models W16 and WJC-2 have quite different scalar off- 
shell coupling constants, but yield the same triton bind- 



TABLE III: Root-mean-square charge and magnetic radii in 
fm. The experimental values are from Ref. [Til ]. 



Model 


r c h( 3 H) 


r ch ( 3 He) 


rw( 3 H) 


r mag ( 3 He) 


W16 (CIA) 


1.718 


1.900 


1.915 


2.001 


W16 (CIA-0) 


1.720 


1.901 


1.915 


2.001 


WJC-1 (CIA-0) 


1.700 


1.879 


1.901 


2.035 


WJC-2 (CIA-0) 


1.722 


1.904 


1.904 


2.027 


Experiment 


1.755 


1.959 


1.840 


1.965 




±0.086 


±0.030 


±0.181 


±0.153 



ing energy. The close similarity of the 3iV form factors, 
at least up to intermediate values of Q, implies that the 
electromagnetic structure of the 3JV bound state is not 
modified too much by the scalar off-shell coupling. This 
conclusion receives even stronger support from the ob- 
servation that also the IARC calculation, which of course 
has no off-shell couplings at all, essentially coincides with 
both W16 and WJC-2 up to about Q = 6 fin" 1 in the 
charge form factors, and up to somewhat smaller values 
of Q for the magnetic form factors. 

It follows then that most of the different behavior of the 
WJC-1 form factors cannot be attributed to the models 
larger scalar off-shell couplings [f| , unless the dependence 
turns out to be highly non-linear. 

What about the contributions of Z-diagrams and pion 
exchange currents? In Ref. |4] it was shown that pion 
exchange currents bring the calculations closer to the ex- 
perimental data. The form of these exchange currents 
depends on the structure of the irNN coupling. For 
pseudoscalar coupling there is no jttNN contact interac- 
tion, but there are large contributions from Z-diagrams. 
The opposite is true for pseudovector coupling. Here, 
the minimal substitution q^ — > q^ — eA into the momen- 
tum dependent interaction, g^NNl^^^m (where q v is 
the pion momentum), leads to the contact interaction, 
e5-7rAfAr7 5 7 A1 /2m, but the Z-diagrams that are produced 
by a pseudovector interaction are strongly suppressed 
(and vanish in the nonrelativistic limit). Furthermore, 
there is an equivalence theorem that has been know for 
many decades 21]: in the nonrelativistic limit, the Z- 
diagrams derived from pseudoscalar coupling are identi- 
cal to the contact interaction derived from pseudovector 
coupling (and the pseudovector Z-diagrams vanish) . The 
pseudovector contact interaction equals the pseudoscalar 
Z-diagram. 

This discussion is helpful in interpreting the difference 
between our results for WJC-1 and WJC-2/W16/IARC. 
The CIA calculations reported here include Z-diagrams 
(to all orders) but do not include contact interactions. 
These are included in the diagrams G-J of Fig. [TJ and 
are excluded from both the CIA and CIA-0 calculations. 
They must be added separately, just as in the work of 
Ref. [H]. Now, the pion-nucleon vertex can be written 
in the general form g^NN^^f 5 + (1 — tyl 5 $ /2m], where 
A is the pseudoscalar-pseudovector mixing parameter [|| . 
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FIG. 3: (Color online) Magnetic form factors of the 3iV bound states, 3 H (first row), 3 He (second row), and the isoscalar (third 
row) and isovector (fourth row) combinations. In each case, the figure on the left shows the form factor in the traditional 
semi-log plot, while the figure on the right shows the same form factor divided by a scaling function || on a linear scale. The 
meaning of the various curves is the same as in Fig. [2] 
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The NN models used in Ref. Q such as W16, as well 
as WJC-2, use pure pseudovector coupling, with A = 0. 
Hence the Z-diagram contributions of these models are 
very small, and it is not surprising that they are quite 
close to the IARC result. However, in WJC-1 the neutral 
and charged pions are treated separately, and the mixing 
parameter for the charged pions is X 7I ± = —0.312. The Z- 
diagram contributions from this model should be large, 
but of the opposite sign from those from a pure pseu- 
doscalar theory (corresponding to X n ± = +1). Hence 
we can expect the Z-diagram contributions from WJC-1 
to move the theory further away from the data, which 
is what we observe. We expect this effect to be more 
than cancelled once the contact interactions are included, 
which for WJC-1 will have a strength 1.312 times a pure 
pseudovector coupling. 



It is certainly not possible to draw very strong conclu- 
sions based on these results alone. There are simply too 
many variables in play, and it would require many more 
test calculations to try to disentangle them. However, it 
is a very interesting situation to have essentially on-shell 
equivalent interactions, which even agree in the 3iV bind- 
ing energy, but lead to different 37V form factors. It has 
often been argued that electromagnetic probes provide 
a means to distinguish otherwise equivalent NN inter- 
action models, and in this case we can actually see it 
happening. From this point of view, it is perhaps less 
surprising that WJC-1 differs somewhat from the other 
models, but rather that WJC-2, W16, and IARC are so 
close to each other. After all, IARC is calculated from the 
largely phenomenological nonrelativistic Argonne AV18 
two-nucleon and UIX irreducible 3iV force, where the 
latter is needed to make up for the missing 3N binding 
energy of the AV18 potential alone. In contrast, the CST 
models do not add any irreducible 3iV forces, and it is 
through off-shell couplings — purely relativistic effects — 
that effective 3.ZV forces are implicitly generated. More- 
over, relativity is implemented in very different ways. It 
is not at all obvious that the two approaches should yield 
so similar results. 



It would be premature to favor one or the other of the 
models WJC-1 or WJC-2 at this time. Their ability to 
reproduce the data can only be judged rigorously after 
all — or at least the dominant — interaction currents are 
added. Compared to the interaction currents of Ref. [J], 
some are already accounted for in CIA, others need to be 
added as well, such as the boson-in-flight terms. Prob- 
ably more important than the latter, there are interac- 
tion currents induced by total-momentum dependencies 
in the vertices due to off-shell boson-nucleon coupling, 
which have never been evaluated. It may very well turn 
out that they have a stronger effect with WJC-1 than 
with WJC-2. 



IV. CONCLUSIONS 

We have performed the first calculations of the electro- 
magnetic 3iV form factors with the new covariant two- 
nucleon interaction models WJC-1 and WJC-2, which 
yield an excellent description of the neutron-proton ob- 
servables below 350 MeV for the most recent 2007 data 
base @. The form factors were calculated in Complete 
Impulse Approximation (CIA), in which — for practical 
reasons — we replaced 3iV vertex functions with two off- 
mass-shell nucleons by corresponding vertex functions 
with only one nucleon off mass shell. This procedure 
of approximating the full CIA results, denoted as "CIA- 
0", was tested with the older two-nucleon model W16, 
for which the full CIA result is also available, and found 
to be of very good quality. 

We compare the form factors of WJC-1 and WJC-2 to 
those of W16, and also to calculations of nonrelativis- 
tic impulse approximation with relativistic corrections 
(IARC) by Marcucci and collaborators [3, [nj. Relat- 
ing the observed differences in the various results to the 
underlying nuclear dynamics, we reach the following prin- 
cipal conclusions: 

(i) The 37V binding energy determines the electro- 
magnetic 3-ZV form factors in impulse approximation up 
to unexpectedly large values of the transferred momen- 
tum. The closely related family of models investigated 
in Ref. |2j , where variations in the strengths of the scalar 
off-shell coupling led to significantly different 3iV binding 
energies, showed much larger changes in the form factors 
than the models considered here, which all have the same 
binding energy. 

(ii) The scalar off-shell coupling does not directly ex- 
ert strong influence on the shape of the form factors. 
When the 3-ZV binding energy is constrained to be equal, 
different scalar off-shell coupling strength can yield very 
similar form factors, as one can see comparing WJC-2 
and W16. 

(iii) In some cases, model WJC-1 deviates moder- 
ately from the others. This appears to be due to its 
mixed pseudoscalar-pseudovector pion-nucleon coupling 
(WJC-2 and W16 have pure pseudovector pion-nucleon 
coupling, whereas in the nonrelativistic framework of 
IARC the two couplings are equivalent). In CST, pseu- 
doscalar pion-nucleon coupling automatically includes Z- 
diagrams, while they are suppressed for pseudovector 
coupling. When Z-diagrams are effectively added to 
IARC in the form of -fnNN contact interactions Q, the 
calculated form factors move closer to the experimental 
data, whereas the WJC-1 form factors lie further away 
than the models with pure pseudovector coupling. This is 
consistent, because the sign of the pseudoscalar coupling 
in WJC-1 is opposite to the one used in Ref. [J]. 

(iv) The results of this work confirm the conjecture for- 
mulated in the Introduction, namely that the reason for 
the good agreement of the CST models with the IARC 
results is the the suppression of Z-diagrams through pseu- 
dovector pion-nucleon coupling. 
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(v) The CST two-nucleon interaction models WJC-1 
and WJC-2 not only give an excellent fit to the avail- 
able two-nucleon scattering observables, but also provide 
a solid basis for a relativistic theory of the 3N system. 
Without additional irreducible 3N forces, the 3N binding 
energy is reproduced, and the electromagnetic 37V form 
factors turn out very similar to previous nonrelativistic 
results. No unusually large interaction currents seem to 
be required in order to achieve a quantitative description 
of the experimental data. 
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